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We consider the motion of an axisymmetric column of Navier-Stokes fluid with a free surface. Due 
to surface tension, the thickness of the fluid neck goes to zero in finite time. After the singularity, 
the fluid consists of two halves, which constitute a unique continuation of the Navier-Stokes equation 
through the singular point. We calculate the asymptotic solutions of the Navier-Stokes equation, 
both before and after the singularity. The solutions have scaling form, characterized by universal 
exponents as well as universal scaling functions, which we compute without adjustable parameters. 
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I. INTRODUCTION 



The breakup of free-surface flows has been an object of intense research from the advent of hydrodynamic theory, 
and in particular the discovery of surface tension Namely, surface tension is the driving force behind this 

phenomenon, as it tends to reduce the surface area by decreasing the radius of a column of fluid. This indeed leads 
to the formation of drops, as is seen most clearly from Rayleigh's |^ stability analysis of an infinite cylinder of fluid 
with radius tq. 

He considered perturbations of different wavelengths and calculated their growth rates. While long wavelength 
perturbations result in the smallest surface area, they require large mass transport between maxima and minima. 
Both effects strike a balance at the wavelength A « 9ro, corresponding to the fastest growing mode. This type of 
analysis subsequently has been greatly refined, for example including viscosity surface charges |6|, or higher order 
nonlinear effects [Q. 

However, even higher order perturbation theory rapidly becomes inadequate as the thickness of the fiuid neck goes to 
zero at a point, and fluid is expelled from this region with increasingly high speed. Near the singularity, characterized 
by a blow-up of local curvature, and of the velocity at the pinch-point, nonlinear effects will soon dominate the 
dynamics. An asymptotic scaling theory of this singularity, where surface tension, viscous, and inertial forces are 
balanced, has been presented very recently jsj. 

But eventually the size of the neck or the times scale on which it is moving will reach microscopic scales, and a 
hydrodynamic description breaks down altogether. For example, the neck will evaporate somewhere close to the pinch 
point, where it has minimum thickness. Shortly after that, new surfaces will have formed on either side, and this time 
the dynamics is described by two separate Navier-Stokes problems. The physical question we address here is whether 
the new initial conditions depend on the microscopic mechanisms behind the breakup. In other words, taking two 
different kinds of fluids with the same surface tension, density, and viscosity, will the breakup look the same on scales 
larger than the microscopic ones? 

We will indeed show that drop formation is a hydrodynamic phenomenon in the above sense. Namely, we construct 
asymptotic solutions to the Navier-Stokes equation after breakup, which describe two separate surfaces and which 
are unique continuations of the solutions before breakup. The physical origin of this uniqueness lies in the properties 
of the solution before breakup |^ . The diameter of the fluid neck does not go to zero uniformly, but only inside a 
"hot" region around the pinch point. Outside, the solution is static on the time scale of the central region. As one 
approaches the singularity, the size of the hot region goes to zero. Hence by the time microscopic mechanisms become 
important, their action is confined to an extremely small region in space. The continuation is achieved by matching 
the outer parts of the solution before breakup onto the corresponding regions after breakup. Since the outer parts 
are virtually unaffected by the microscopic dynamics, this procedure yields universal continuations. 

This seems to be the first example of a partial differential equation uniquely describing a "topological transition" 
1^. The result is also important for numerical simulations, which usually rely on some ad-hoc prescription for the 
formation of a new surface pO| , pT[ | , or for breakup in related physical situations . 

Our paper is organized as follows: In Section 2 we derive a one-dimensional approximation of the Navier-Stokes 
equation | , valid as the ratio e of the radial to the longitudinal scale of the flow is small. They have self- similar 
pinching solutions, which are described by a pair of scaling functions (/){£,) and tp{£,) for the radius of the fluid neck 
and the velocity, respectively. As the time distance from the singularity goes to zero, the slenderness parameter e for 
this solution vanishes, making it an exact solution asymptotically. 

The scaling functions (j) and tp obey two pairs of ordinary differential equations, one for the time before breakup, the 
other for the time after breakup. For most of the rest of this paper, we will be constructing unique solutions to those 
equations. In the third section we consider the similarity equations before breakup. Shortly before the singularity, the 
fluid far outside the pinch region is no longer able to follow the motion near the pinch point. This leads to boundary 
conditions for the similarity functions at inflnity, and together with a regularity condition in the interior, a unique 
solution of the equations is selected. We compute this solution numerically. 

The same procedure is adopted in the fourth section for the similarity equations after breakup. Here the solution 
is matched onto the profiles before breakup. This solution has two halves, each of which is fixed uniquely by the 
matching. 

The concluding discussion gives an example for the breakup of a real fluid, which could be measured experimentally. 
We also supply numerical evidence for the uniqueness and stability of our theoretical predictions, and discuss related 
work. 
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II. SIMILARITY EQUATIONS 



Let us begin by formulating the Navier-Stokes problem for an axisymmetric column of fluid, where we assume the 
azimuthal velocity to be zero. A sketch of the geometry of the problem can be found in Figure 1. For a fluid with 
kinematic viscosity v, surface tension 7, and density p the Navier-Stokes equation reads in cylindrical coordinates 




dtVr + VrdrVr + VzdzVr = —drp/p+ l'{d'^Vr + S^Wr + BrVr j T — VrjT^)-, (1) 



dtVz + VrdrVz + VzdzVz = ^B^vIp + viplvz + dlvz + drVz/r) - g, (2) 
with the continuity equation 

drVr + dzVz + iw/r = 0. (3) 

The acceleration of gravity points in negative z-direction. Here Vz is the velocity along the axis, Vr the velocity in the 
radial direction, and p the pressure. There are two boundary conditions, coming from the balance of normal forces, 

nan = -7(l/i?i + l/i?2), (4) 

and tangential forces 

n (7 t = 0. (5) 

In we denoted the outward normal and tangent vector to the surface by n and t, a is the stress tensor, and 

(l/i?i + l/i?2)/2 the mean curvature. A standard formula for bodies of revolution gives 

1 1 1 dlH 

(6) 



i?i i?2 H{i + {dzHYY/-^ (l + (a,i?)2)3/2' 
where (z, t) is the radius of the fluid neck, as seen in Figure 1. The equation of motion for H{z, t) is 

dtH + VzdzH^Vr\r=H, (7) 

which says that the surface moves with the fluid at the boundary. 

Equations constitute a complex moving boundary value problem, which we want to investigate near a 

singularity, where nonlinear effects are bound to become dominant. The reason exact solutions, valid arbitrarily close 
to the singularity, can nevertheless be found, is that only very few terms in the equations contribute to the leading 
order force balance. Thus to proceed, we first have to identify those leading order terms. We will then construct 
explicit solutions to the leading order equations and demonstrate their consistency with both the internal structure 
of the Navier- Stokes equation and with boundary conditions. 

The relevant terms are identified using two properties of the singularity to be validated later: 

(i) The singularity is line-like, i. e. its axial extension is much greater than its radial extension. 

(ii) Surface tension, viscous, and inertial forces are equally important near the singularity. 

Conditions (i) and (ii) are now incorporated into a perturbation theory. According to (i) we will assume that the 
motion of the fluid at a given time is characterized by an axial length scale Iz and a radial length scale for which 

tr = ttz. (8) 

where e is some small parameter. The physical meaning of e will come out later from the description of the singularity. 
Also introducing a time scale tz of the singularity, we can nondimensionalize all quantities according to 

r — irf , z — £zZ , t — tzt , (9) 

iz P tip 

tz p Hp 

t o 5 ' ^ t2 ^ ' 

Lz p Lz P ^z 
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The scales iz, ir, and are defined to be constants, so their derivative with respect to time is zero. However, 
one must bear in mind that the characteristic scales of the singularity change, so £z, ir, and tz will be different in 
different stages of the singularity formation. Since there are two length scales £z and ir, there is a certain freedom 
in the nondimensionalization of the material parameters v, j/p, and g. This freedom is completely specified by the 
exponents n, m, and I in (^). We will see below that the exponents are fixed by the requirement (ii). 

Since the radial extension of the fluid is small, we can expand all fields in the dimcnsionless radial variable r : 

oo 

= ^S2j(i,i)(ef)2j' , (10) 

1=0 



Vr{z,f,i) 



2j + 2 



(11) 



and 



p{z, r, f) = ^ p2j {z, i) {erf^ 



(12) 



The definition of Vr automatically ensures incompressibility. We now insert (H)-(^2]) into the equations of motion 
(|l|)-(0), and compare powers in e. The lowest order expressions result in a closed set of equations for vq and H, 



H 
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(13) 



diH + vod-zH = -idsvo)H/2 . (14) 

To obtain closure at higher orders in e, one needs to expand each of the coefficients V2j and ■p2j, as well as H into a 
separate power series in e. There then exists a consistent representation of (|l|)-(0) to all orders in e ||l6|. We will not 
be concerned with the explicit form of the higher order equations here, so for simplicity we use the notation vq and 
H (or its nondimcnsional counterpart) for the lowest order terms in the expansion in e. 

It is evident from ( p^ that the exponents to, n, and I determine the balance of forces at leading order. Since the 
l/H term, which comes from the radius of curvature perpendicular to the axis, is driving the instability, it must 
clearly be present and in fact becomes infinite at the singularity. At the small scales involved in singularity formation, 
viscosity will also be important. Finally, velocities are expected to blow up as ever smaller amounts of liquid are 
driven by increasingly large pressure gradients. Hence we also expect inertial effects to be involved asymptotically. 
Since the acceleration of the fluid diverges at the pinch point, the constant acceleration of gravity will drop out of 
the problem. This is precisely the assumption (ii), incorporated by choosing to = 1, n = 0, and / > in (^, which 
leads to an equation where surface tension, viscous, and inertial forces are balanced, while gravity is irrelevant. These 
assumptions will be tested for consistency later. 

We now identify the scales involved in the formation of the singularity. It is crucial to notice that all external 
length and time scales, which are imposed by boundary and initial conditions, do not enter the description of the 
singularity. In a jet experiment, for example, external scales would be the radius of the nozzle and the period of the 
driving frequency. 

Near the singularity, the length scales characterizing the solution become arbitrarily small, while time scales become 
shorter and shorter as one approaches the singularity. Hence the singularity moves on scales widely separated from 
the external scales. It is for this reason that for the mathematical analysis of the singularity we do not have to make 
the boundary or initial conditions explicit. Boundary and initial conditions will become important when we describe 
numerical simulations of real experiments, which confirm the consistency of our approach. 

The proper units in which to represent the motion near the singularity can thus involve only internal parameters 
of the fluid. This leaves us with the units of length and time 

t. = {pi^^)h , = (p'^^^)/7' • (15) 

Assuming that the singularity occurs at a point zq, and at a time to, the space and time distance from the singularity 
is properly measured as 
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z' = (z-zo)/^. , (16) 
t' ^{t- to)/U ■ 

The units £1, and are a measure of the width of the critical region, and are fixed for a given fluid. Singular behavior 
is expected for \z'\ <C 1 and \t'\ <C 1. Note the conceptual difference to the characteristic scales £z, £r, and tz of the 
singularity, which change in time. 

In the variables z' and t', the fluid velocity and the neck radius are: 

«(z',t')^^«o(^,t) , (17) 

h{z',t') = e-^H{z,t) . 

Keeping the same terms as in (|l^),(|l^ with m = 1, n = 0, and / > 0, we find in the limit e ^ 

dt'V + vdz>v = -dz' l-j +3 ^—^ , (18) 

dt'h + vdz,h^ -{dz'v)h/2 . (19) 

All material parameters have dropped out of the equations, since everything has been expressed in units of £1, and t^. 

At this point it is worthwhile to pause and notice that we have already succeeded in reducing the original Navier- 
Stokes problem in two spatial dimensions and in time with a moving boundary to just a coupled set of equations in 
one space dimension and time, at least for small e. Approximations for thin liquid threads of the type described here 
have in fact a long history, see ||l^ for a (by no means complete) list of earlier references. However it seems that 
(|l8|) , (|l9|) , which contain the correct surface tension, incrtial, and viscous terms, were first derived in ||l^. Another 
related approach goes by the name of Cosserat equations, see for example pT) . In all previous work except |^ though, 
the resulting one- dimensional equations are treated as model equations, whose quality of approximation depends on 
the particular physical situation for which they are used. In the present paper, we will show that (|l^), ([l^) become 
exact close to pinch-off. 

To this end we have to identify the parameter e. From the definitions (llq) we find 
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tz/U^e^^ , e./e.^e^ . (20) 

Thus up to constants is the characteristic time scale of the singularity, written in units of t^. But the only such 
time scale is the nondimensional time distance from the singularity \t'\ itself. Hence \t'\ serves as the desired smallness 
parameter. We introduced the modulus of t' here, since we need a measure of the time distance before and after the 
singularity. As \t'\ — s- 0, all higher order terms vanish and only the leading order equations (18),(p^) remain. By the 



same token, the axial and radial length scales behave like £z ^ £u\t'\^^'^ and £r ~ £v\t'\- Thus close to the singularity, 
all length scales become arbitrarily small compared to any external length scale, just as we asserted above. 

As a corollary to this absence of any fixed length scale in the problem, we expect singular solutions to have the 
similarity form 

h^\t'rm , (21) 
v = \t'rm , 



where the similarity variable is defined as ^ = A similar ansatz has been used in |jT^ for a study of inviscid 

flow, but in a different geometry. 

The values of the exponents ai, 02, and (3 are inferred immediately from dimensional analysis. Namely £r ^ £u\t'\ 
implies ai = 1, £z/tz ~ i£u/t„)\t'\~^^^ is a typical velocity scale, giving 02 = —1/2, and /3 = 1/2 follows from 
£z ^ £v\t'\^/'^ . The appearance of fractional powers forces us to use the modulus of t' in the scaling laws (^l]). The 
type of similarity solutions we are going to investigate is thus 
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h=\t'\m , (22) 

v = ±\t'\-^/'m , 

C = ±z'/\t'\'^' . 



The two different signs take care of identical solutions with different parity. The acceleration of the fluid diverges like 
ji'l^"^/^, and surface tension, viscous, and inertial forces are balanced. Since \t'\ ~ we conclude that the exponent I 
in (|9|) is Z = 3, which is consistent with our previous assumptions. 

Inserting (|2^) into ( |T8| ) and (^9|) we find that the asymptotic equations of motion indeed have scaling solutions, 
where the scaling functions (j) and -0 obey the equations 

s(V'/2 + CV72) + i'lp' = (t>'l4>'^ + aV'" + 4>' I (j) , (23) 



s{-4> + ^072) + ■4'4>' = ~V''0/2 . (24) 

The prime refers to differentiation with respect to ^. The terms in brackets come from the time derivative, s = 1 
refers to the time before the singularity (i < to), s = — 1 to the time after the singularity {t > to). 

Hence close to the singularity, \t'\ <C 1 and \z'\ <C 1, we have further reduced the problem to a set of two ordinary 
differential equations. To find unique solutions of ( p^) and ( p^ we still need to formulate appropriate boundary 
conditions. This and the numerical integration of (|23|)7(p4[) will be the subject of the next two sections, first for t < tp, 
and then for t > t^. 



III. BEFORE BREAKUP 

In this section we consider the similarity equations (^), ( |2^ ) for s = 1, i.e. before breakup. Some of the calculations 
relevant for the next section will be done for general s. We show that the similarity equations have precisely one 
physically allowed solution, and compute it. Therefore singular solutions are completely universal: once the origins 
of the space and time axes are fixed by specifying zq £ind to, there arc no more free parameters. The relevant units of 
length and time are set by the fluid parameters. 

As the similarity equations are of first order in 4> and of second order in -0, solutions are specified by three initial 
conditions 4i{£,i), '0(CO, i^'i^i) at a reference point ^i. Universality implies that we need to find three conditions 
which uniquely fix the physically allowed solution. 

For the first condition, suppose we choose a small region of width £i,d around the singularity, such that S <^ 
mm{l,L/£^), where L characterizes some outer length scale. For \z'\ < 6 and \t'\ <C 1 we are well within the critical 
region of the singularity, and effects of the boundaries are negligible. Thus the similarity equations (p3|),(p4[) apply for 
\z'\ « S and we have |t'|(/()(±5|i'|~^/^) « h(±S,t'). First we observe that the point \z'\ = S goes to infinity in similarity 
variables as \t'\ 0. Second, in this limit h at \z'\ = S will not be able to follow the motion of the singularity, whose 
width decreases like |t'|^^^, and whose time scale goes to zero with \t'\. Hence h{±6,t') must approach a finite value 
as \t'\ —^ 0. To be consistent with this physical requirement, </)(^) must grow quadratically as |^| goes to infinity. 



Hence we have two conditions on the solutions of (23),(|24|): 

a-) "0(0 need to be regular on the real axis C ~ oo, +oo[. 

b) For ^ — !■ ±oo, (/>(^)/^^ should approach a finite limiting value. 

Conditions similar to b) have also been employed in |p^ . Note that the physical concept behind our argument is 
inertia, which prohibits the large amount of fluid far away from the singularity to move with the fluid in the skinny 
pinch region. 

We will now show that the requirements a) and b) completely determine the solution of the similarity equations. 
In particular, we do not have to specify the limiting values of 0(0/^^, they rather come out of the solution of the 
problem. This is consistent because in our analysis we deal exclusively with the equations of motion valid close to the 
singularity. No input from regions where the expansion is not valid is needed. Thus boundary or initial conditions 
can enter the problem only implicitly, as they determine the position of the singularity zo,to. 

Let us begin by examining the behavior of solutions for ^ ±oo. It is advantageous to first eliminate from the 
problem, leaving us with a third-order equation for i/j po|. To this end (p3) is written as 
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s - V>72 

On one hand this equation can be used to express (j) in terms of tp, 

3V^")/i^ - 6^']-' , (26) 

where we have introduced the notation 

i^ = (s-V72)/(^/' + sC/2) , (27) 
/ = (sC^ + V'')/2 . 

On the other hand, writing ip as an integral over the kernel K we have 

= 0(eo)exp|^^X(C)dc| . 
Inserting this into (|2^), taking the logarithm, and differentiating, we find 

= ^ { (I'/K)' + I' + {K'/K^ - 3) - &^'K] , (28) 

which is a single equation just in terms of ip. 

Plugging the ansatz -0 = B^" into (p8|), one finds the leading order behavior on the right hand side to be 



2B 

Thus i/" must decay like I/^ or 1/^"^ at infinity for the terms to cancel. In particular, growth of ijj at infinity is 
prohibited, since the "inertial" term /' is quadratic in ip, and would grow faster than any other term in (p8[). 
Thus one is lead to an asymptotic expansion of the form 

oo 
^ 1=0 

Only odd powers appear, since (^8|) is invariant under the transformation ^ — ^ and ip —>■ —ip. Using ( p6| ) we can 
calculate the leading behavior of (j) corresponding to ([29|): 

= aoe'[l + O(r')] , (30) 
ao = 2/[66o - s6i - 6o] ■ 



This means ( P9|) represents precisely the physically relevant solutions we are interested in. To further investigate 
the expansion (E9), we derive recursion relations for the coefficients bi to arbitrarily high order. The lowest order 
expressions are 

b2---^[3bl + 7sboh] , (31) 

63 = f [-ml + 9b^ ~ USsbobi - 9sblbi - lObj - 862 (10 + 36o)] ■ 
8 

All bi are thus determined by just two free coefficients, 60 and 61, or by virtue of (^), qq and bg. However, the 
expansion ( |2^ ) is only asymptotic, as for large i the bi grow like 

- {-12)H\ . 

This means for large ^ all solutions of ( |2^ ) are up to exponentially small correctionssiyen by a two- parameter 
family of functions ipagboiO (111' '^hich behave like 1/^ asymptotically. The expansion ( p9| ) is asymptotic to ipaoboiO 
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and for sufficiently large ^ can be used to compute 4'aobo (0 ^o any desired accuracy. To understand the significance 
of this observation, we have to investigate the stability of the functions V'aofco(0- 

Doing so turns out to be slightly more convenient in the original space of initial conditions ''/'(Oi V^'(C)); 

where ^ 3> 1 is kept fixed. Denoting by (j)ab{C) the function (j) corresponding to ijjabiC)^ s-'"*' interested in particular 
in perturbations which carry us out of the two-dimensional manifold of initial conditions {4>ab{S,),'4^ab{C)i''P'abi^))- 
Differentiating with respect to a and 5, we find that to leading order in ^, (0, 0, 1) is a vector normal to this manifold. 

We now consider small perturbations relative to the solutions 4'ab,'4'ab'- 

m^4>ab{m + ^m) (32) 
m = ^abmi + ^2{0) 
^\i) = ^'ab{m + ^m)- 



The correction e3(^) describes the behavior of perturbations perpendicular to the plane of asymptotic solutions ~ 
-0 ^ Inserting (|2|) into (^), and linearizing in the ti reveals that to leading order £3 behaves like 

e3(e) = e3(aexp||(e-e)} . (33) 

Hence for s = 1 an arbitrarily small perturbation introduced at ^ will carry the solution away from the physically 
relevant manifold as |^| oo. Only a two-dimensional manifold of solutions is consistent with ~^ const as ^ 

tends to -|-c» or — oo. This means the requirement b) corresponds to two constraints on physically relevant solutions. 
Since the equations are of third order, we need to find one additional constraint to uniquely fix the allowed solutions. 
It is worth remarking that the unstable growth (|3^ ) comes from the presence of the viscous term ip" in (|2^). Hence 
in a strictly inviscid theory no selection would take place. 

To find the third constraint we look at condition a), saying that cj), ip be regular. Considering ( p5|) this is a nontrivial 
condition, as tp must be bounded and hence there is a point with 

V'(Co)+eo/2 = . (34) 

Therefore, since s = 1, the denominator in ( ]25| ) will vanish at leading to a singularity unless the condition 

^'{^o) = 2 (35) 

is also met. To explore the corresponding regular solutions, we expand i/j around ^q- 

oo 

m^Y.'^^^^^^oT ■ (36) 



The function (j) can again be recovered from (E6|). We find 



do = -Co/2 , (37) 
di = 2 , 

5^0 



8(3-l/0o) 

(104- 16560o)rfi/75-H60o 
2 - 36(j)o 
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where the first two equations follow from ( |34D and (^5|). Just as in the expansion around ^ = ±oo, all coefficients di 
are determined by only two coefficients, (t>o — 0(Co)- We verified this statement by deriving recursion relations 

for the di to arbitrarily high order. This time the expansion has a finite radius of convergence, whose value depends 
on the initial conditions S,o,4'o- 

It is worthwhile to comment on the physical significance of ^o- The equation of motion for the position z'^ of a 
marker on the surface h is 

dtz',{t) = v{zi{t),t) . (38) 
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Rewriting z'^ in similarity variables, = \t'\ ^I'^z'^^ and measuring time on a logarithmic scale, s = — In we find 

a,e.(s)-6/2 + V^(6) , (39) 

which is the convection equation in similarity variables. Hence at the point a-s defined by (|3^), a surface marker 
on is at rest. Regularity properties on such "stagnation" or "sonic" [g2[ points often play a similar role in selection. 

To explicitly compute the unique solution of the similarity equations before breakup, consistent with a) and b), 
we proceed as follows: we choose a pair (^Oi'/'o) a-nd compute the Taylor coefficients di to sufficiently high order. 
This leaves us with a series representation of V' in a disk around ^o- From there onwards, ( psj ) has to be integrated 
numerically. Since as |^| — > oo solutions must be exponentially close to a two-parameter family of functions il)ah 
which are "repellent" solutions will generically not extend to infinity, but rather end up in a singularity at finite ^. 
Dominant balance in ( |2^ ) reveals that those singularities have the leading behavior ~ — f)""'^. Only a 
one-dimensional submanifold in (^oi<?^o) is consistent with the solution extending to either -|-oo or — cx3. The point 
(foj 00 ) where both cross corresponds to the unique solution we are interested in. 

Our numerical procedure was to introduce and ^~ as the values of |^| where exceeded a certain bound as 

^ — > cx) or ^ ^ — cx), respectively. We then optimized and_(/)o to give maximum values of and As solutions 
deviate exponentially from '^ab^ the "window" around (^o,0o), which allows for solutions extending up to a given 
1^1 gets small very rapidly with |^|. Thus this method allows for a very accurate determination of and 00 . The 
numerical values we found are quoted, together with other characteristics of the solution, in Table 1. These results, 
with the inclusion of the asymptotic expansion (p9|), now allows us to plot the scaling functions before breakup, c/)"^ 
and V''^, in Figure 2. 

As seen in Table 1, the stagnation point is extremely close to the point ^mm where (\)^ is minimum. This means 
that in the frame of reference of the surface, fluid is expelled on either side of the minimum. From z^j„ = |i'|^^^Cmm 
one sees that the minimum moves with velocity Wmm = (Cmm/2)|t'|~"'^^^. 

To make contact with the qualitative description of the singularity given in the Introduction, we schematically 
divide the similarity solutions into three regions: A central region around the minimum of size £,centrah say, where 
(j) is almost constant, and outer regions on either side, where </> is quadratic. In this simplified picture, in physical 
space there is a region of size £,centrai\t'\^^'^ around zq, where the diameter of the neck decreases linearly in time, and 
the velocity diverges like Outside this region, both the thickness of the fluid neck and the velocity field are 

constant. Hence as \t'\ — > 0, at any given point z 7^ zq the solution will become static, and the singularity only occurs 
at a point zq in space. In terms of some microscopic length imicro, one can estimate (molecular) mechanisms to be 
important in a region of size £,centrai{^micro('vY^'^ ■ 

But perhaps the most striking feature is the extreme asymmetry of f/)"*" and ij)^ . Indeed, the values of Oq , describing 
the amplitude of </)"*" as ^ — > ±00, differ by almost four orders of magnitude. Intuitively, an asymmetric solution is to 
be expected |Q. Namely, pressure will be higher in the slender part of the solution, pushing fluid over to the right. 
This will cause the right side of the solution to fill up with even more fiuid and get steeper. Eventually, this mechanism 
is only checked by viscosity. But this argument does not even give an order-of-magnitude estimate of /a^ . So 
clearly there is the need for a fully analytical theory of the selection problem, which gives at least reasonable estimates 
for the numbers in Table 1. 

Another, perhaps related problem pertains to the uniqueness of the above solution. In principle, the one-dimensional 
submanifolds corresponding to the correct asymptotic behavior as ^ — > -l-oo and ^ — > —00 could have several crossings, 
giving a discrete family of solutions. The most reasonable guess for a different form of solution would be a symmetric 
one, which would then be highly unstable, since small asymmetries would amplify according to the above mechanism. 
Since Co = for such a solution, 0o would be the only free parameter, which needs to be consistent with the behavior 
at infinity. We carefully looked for solutions of this type, but found none. Therefore, to the best of our knowledge, 
there is precisely one possible solution, but for a final word we must await a rigorous mathematical theory. 



IV. AFTER BREAKUP 



We now turn to times t > to, i.e. after breakup. In terms of the similarity equations (|23|),(24) this means we have 
to put s = — 1. But apart from the difference in the equations, there is a completely new type of problem occurring 
now, related to the mathematical description of a receding tip. 



To understand this, let us consider the asymptotic equations (18),(|19|), which contain the leading order terms of the 
Navier-Stokes equation as the slenderness parameter e goes to zero. But this description breaks down as one reaches 
the tip, which is assumed to be at z[^^{t), see Figure 3. Namely, the slenderness assumption means that dz'h is of 
order e, while dz'h actually diverges as z' z'^^^. Indeed, both the asymptotic form of the pressure gradient {dz'h)/h? 
and of the viscous term dzi[{dziv)h?]/h'^ diverge as /i ^ and dz'h — > 00. 
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On the other hand, the complete Navier-Stokes problem has no singularities as long as \t'\ > 0. Surface tension will 
ensure that the gradient of the curvature remains finite. Hence there is a small region around the tip, whose width 
goes to zero as e ^ 0, where higher order terms in the Navier-Stokes equation will be important. Its size iup can be 
estimated by saying that the asymptotic equations become valid as dz'h becomes of order unity at the edge of this 

(iiy\t'\, using the known scaling of the radial length scale £r with \t'\. 



region. Thus, since dz'h « Ir/^upi we have tup 

Now we transform to similarity variables ^ = z'/jf'j^/^, where £,tip = -Ztip/I^T^^ is the position of the tip. Since the 
width of the tip region shrinks as \t'\, it will go to zero like even in ^-variables. In the neighborhood of any 

C G]6ipj oo[ the similarity equations will be valid as \t'\ 0. Thus to capture the leading self-similar behavior of the 
Navier-Stokes equation after breakup, one just has to find the correct boundary conditions for (f) and ij) at S^up- This 
situation is reminiscent of the boundary condition at ^ = ±oo before breakup: For ^ the range of validity of 
the similarity equations extends to infinity, so supplying boundary conditions at ^ = ±oo suffices to uniquely solve 
the problem. 

To derive the correct boundary condition, we proceed as follows: We supplement (^8|),(|l^) with higher order terms 
in e, which regularize the equations at the tip. The corresponding similarity equations in cj) and ijj now still contain 
t' as a parameter, but are finite as ^ — > S^up- This means solutions of those equations can be supplemented with the 
natural boundary condition 4>{S,tip) = 0. 

Then we derive a simplified version of the equations valid at the tip, which we can integrate explicitly, using 
4'i^tip) = as a boundary condition. Now we can take the limit — s- or e — > 0, which leaves us with the correct 
boundary conditions for (p and V'j valid for t' = 0. We also show that this result is independent of the particular 
regularization we have been using, so the result is unique, as expected from the above argument. Once the boundary 
condition has been found, we can solve the similarity equations to find a unique solution after breakup. 

No knowledge of the fiuid motion in the tip region of size £^\t'\ is needed to calculate the self-similar part of the 
solution It remains an interesting open problem to devise a method to compute an approximate solution in the tip 
region. However, since this region becomes arbitrarily small as \t'\ — > 0, we will not be concerned with this question 
in the present paper. 

To construct a regularized version of (Qq) , we observe that it can be generalized in the form p3| : 



dt'V + vdz'V = —dz'P + 



dz' [{dz'v)D' 



(40) 



dt'h + vdz'h^ -{dz'v)h/2 



(41) 



with 



1 

2h 



dE 



dE 



dh dz' d{dz'h) 



Here E = E{h, dz'h) is a surface energy and D — D{h, dz'h) a dissipation kernel. This nomenclature is motivated 
by the fact that dz'P may be written as 



dz'P 



1 d 
dz' 



h'p + idz'h) 



dE 



d{dz'h) 



E 



and hence we have the conservation equation 

d 



dt 
d 



- [h'v^l2 + E{h,dz'h)\ = -{{dz'v)Df 



{v^/2+p)h'^v - v{dz'v)D^ + (dt'h) 



dE 



d{dz'h) 



(42) 



So apart from a surface term this equation says that the sum of kinetic and potential energy decreases with a 
negative definite dissipation rate T) — —{{dz'v)D)^ . In the present context, (^0|),(pl]) are phenomenological equations. 
There are certainly other higher order correction terms present in the Navier-Stokes equation, which have not been 
included. However, the only important point here is that E and D can be chosen such as to make the equations finite 
at the tip. In jlj] we already introduced a variant of (^) , ([il]) with 



E{h,dz'h) ^ 2h{l + (dz'h)^)^/"^ 



(43) 



10 



This energy is proportional to the surface area and arizes naturally when keeping the complete curvature term in the 
boundary condition (jj). 

If the surface at the tip is non-degenerate and the velocity field is regular, we simply have h{z,t) — hQ(t){z' — 
Z(jp)-^/^ + 0(z' — Zjjp)-^/^ and v{z,t) = vo{t) + vi{t){z' — z'^^^) + 0{z' — z'f^p)"^ . As is verified by inspection, the particular 
form (^) of E succeeds in keeping dz'P finite as z' z^^^. Introducing 

D{K d.'h) = /i(3/(l + {d,'hf)Y'^ (44) 

for the dissipation kernel, the same is true for dz'[{dziv)D'^]/h'^ , hence all terms in ( ]40| ) are now finite at the tip. At 
the same time, the asymptotic equations (18),(|T9|) are recovered for e ^ 0, as this corresponds to 

Easymp — (45) 
-^asyrnp — \/3/i. 



By construction, all allowed functions E and D must have the same limit (45). We now insert (|2^) into the 
regularized equations (40),(^). For t > to, denoting \t'\^^'^ by i, we obtain: 



and 



where 



(t> - + ip(t>' = , 



G 



(1+ £2^/2)1/2 £20/2)3/2 ' 



(46) 

(47) 
(48) 



1 + £20'2 



1/2 



Here for simplicity we have used the special forms (^3|) and (^J) for E and D. For £ = we recover the asymptotic 
equations (23),(p4), while for finite I all terms are regular at the tip as 4> and '0 behave like 



To focus on the tip region, we introduce the rescaled fields cj) and ip: 

^(C) = ri[V'(^C + W-6.p/2] , 

In rescaled variables, the equations are 

^2^2 [_^/2 _ ^^'/2 + ^4,'] - ^,^2^,,p/4 = {-G + i^'D^y 

and 

- C^'/2 + = - W2 , 

with 



(49) 



(50) 



G = 



D = 



(1 + ^'2)1/2 (1+0/2)3/2 



(51) 
(52) 

(53) 



1/2 
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In (|5l|)-(53) and ( p4[ ) below, primes refer to differentiation with respect to the rescaled variable C 

The only place where ( still appears is in front of the "inertial" terms on the left hand side of (|5l| ) . This is because 
any fixed region ( e [0, Ci] near the tip shrinks to zero in ^-variables as ^ — > 0. But the fluid at the tip should move 
with the boundary, so it is at rest in the frame of reference of the tip. Indeed, since the left hand side of (^) only 
contains lower order derivatives, the limit £ — > is regular at fixed initial conditions for (/), ip at 0. 

Hence by putting = in ( pl|) we obtain a simplified description of the tip region, which is uniformly valid in 
any fixed interval [0, Ci]- Note that implicitly i is still present by virtue of (|50|). Solutions of the resulting equations 
correspond to a very much blown-up version of the tip. Since the solutions are regular at C = 0, we can employ the 
natural boundary condition (f){0) — 0, and from ( |5^ ) we have G(0) — -D(O) = 0. This means ( pi] ) can be integrated to 
give 

G = i^'D^ . (54) 

We now supply appropriate matching conditions, which express the consistency of ( |52| ) and (|^, valid at the tip, 
with the solutions outside the tip. At fixed ^, 0(^) and ip{£,) are finite in the limit £ ^ 0. Since the tip region gets 
arbitrarily small in ^-variables this is physically reasonable, but has also been checked numerically by integrating 
(46),(|4^). Hence we have to require (j),'tjj to behave like « ki and V'(C)/C ~ '^2 for large ^. In view of the scaling 
(50) this makes them consistent with </>(C))V'(C) finite. Inserting (f>{(^) = ni and ?/'(C) = K2C into (^) and (|^, one 
confirms this ansatz to solve the equations, and finds K2 = —2 and ki — 1/6. So, again considering (po|), the lowest 
order terms of (j) and ip as — £,tip) tends to zero are = 1/6 and ip = ^tip/'^ — 2(^ — ^up). In other words, at ^up we 
have the boundary conditions 

<^(6.p) = l/6 , (55) 



This boundary condition implies that the asymptotic shape of is a step function of height |t'|/6 at the poi nt ztip{t). 

It is important to notice that this result is independent of the particular form of regularization ( 43 ) , ([44|) we have 
been using. For example any other term involving h leads to a term i'^cj) and drops out as € — > 0. Another contribution 
dz' h gives </>' and also does not contribute as we finally set = ki . 

It only remains to formulate boundary conditions for ^ — s- 00. At large distances from the singular point zq both 
the interface and the velocity field should look the same as before breakup. This is the same reasoning that made us 
construct solutions which far away from zq are static on the time scale of the singularity. As the width of the singular 
region shrinks to zero like the large body of fluid outside is not able to follow. Here it provides us with the 

mechanism for unique continuation: For the two solutions to coincide we must require that 



hm 

{— >±oo 

liin V(C)C = 

4— >±oo 



(56) 



We will see that (p5|),(|56D are all the boundary conditions needed to uniquely solve (|23|) , (p4D after breakup. Since 
the constants oq and 60 are different for the left and right hand side of the problem, the solutions will also differ. In 
particular, the value of ^tip consistent with (^6|) depends on oq and bo. The requirement ( ^6| ) thus represents the way 
the properties of the solution before breakup are communicated to the solution after breakup. Inserting the ansatz 



= i/6 + 0i(C-6.p)" + ... , 

^ = 6«p/2 - 2(C - Cup) + eo(^ - Cupf + 



(57) 



into the similarity equations and balancing leading powers we find a = 2/5 and /? = 7/5. We therefore try the general 
expansion 



= Cup/2 + (e - Cup) 



-2 + ^e.(e-6.p)^'+ 



^)/5 



i=0 



(58) 



Again, by (E6h it is sufficient to consider the expansion of ip. The first few coefficients are 
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eo = yf/)! , (59) 

ei =0 , 

120 2 

62 = Tj-(Pl , 

We confirmed, by deriving recursion relations for the to arbitrarily high order, that all coefficients are determined 
by the two free parameters $,up and (pi . Since the power series ( p8| ) has again a finite radius of convergence, all solutions 
starting from £^tip are classified by just two parameters. But for t > to the behavior for |^| ^ cx) is very different from 



the situation before breakup. We now have s = — 1, and according to (33) the asymptotic behavior (/) ~ ip ^ 
is stable. So integrating the similarity equations to infinity, for every value of ^tip and (pi we will find a unique value 
of lim^^oo </>('f)/'C^ and lim^^oo V'(C)C- Hence the boundary conditions ( p6|) are precisely what is needed to uniquely 
fix ^tip and (pi, and thereby uniquely determining the similarity solution (p'~ and ip~ after breakup. 

Obviously, this has to be done for the left and right hand sides separately. The left hand side corresponds to a 
receding neck, the other is the main drop. We denote the values of ^up and (pi by ^„ecfe and (pneck for the left hand 
side, and ^drop and (pdrop for the right hand side. The result of a numerical calculation of and ^~ can be found 
in Figure 4, some of the characteristics of the solution are listed in Table 2. Specifically, the neck recedes with the 
velocity 

S,neck |,/|— 1/2 rcn\ 
Vneck^—^\t\ ' , (60) 

where £,neck/'^ « 8.7. Unfortunately, on the scale of Figure 4 it is hard to see any deviations from a flat interface for 
the drop. Figure 5 below will give a better idea of how the drop is left distorted after breakup. 

It should be appreciated that the unique continuation does not follow from the asymptotic equations (^3|),([2^) 
alone. Rather, we needed to invoke regularity for \t'\ 7^ to derive the boundary condition 4>~{^tip) = 1/6. Indeed, 
( p3| ) and ( [2^ ) with s = — 1 would allow for an infinity of solutions, one for each value of <p{£,tip)- 



V. DISCUSSION 



We have shown that the motion of a Navier-Stokes fluid close to the time of breakup is described by self-similar 
solutions. The corresponding scaling functions, before and after the breakup, are solutions to a set of ordinary 
differential equations. For the solutions to be consistent, both away from the singular point and at the receding tip 
after breakup, boundary conditions have to be imposed. They lead to unique solutions of the similarity equations. This 
means solutions to the Navier-Stokes equation close to the singularity are predicted without adjustable parameters, 
and independent of boundary or initial conditions. It is quite instructive to plot the predicted interface of a real fluid 
at constant time intervals before and after the singularity. Since ti, and are almost on molecular scales for water 
[0, we take a mixture of glycerol and ethanol as a reference fluid, for which (.^ — 72/iTO and ty = 114/is. This is 
large enough for experiments by optical means to be feasible. Measurements of the velocity field are also possible 
p4| . Figure 5 shows three profiles, each 46/is apart, before the singularity (a), and after the singularity (b). This 
corresponds to \t'\ = 1, 0.55, and 0.1. In particular, there is no freedom in the spatial scale of this Figure. The same 
graph should apply regardless of boundary conditions. 

Before breakup, one can clearly distinguish a very slender neck, and the steep front of the adjoining drop. As the 
neck becomes thinner, the minimum moves towards the drop, making the interface even steeper. The greatest relative 
changes in the diameter occur near the minimum, far away the interface is practically static. As one comes closer to 
the singularity, the size of the "active" region, which is still changing, becomes smaller and smaller. 

After breakup, the neck snaps back very rapidly, forming a sharp front at the end. For \t'\ = 1, higher order 
corrections in \t'\ will probably be already important, and the end will look more rounded. As seen in Figure 4, there 
is also some fiuid accumulating at the end in the asymptotic solutions, but this cannot be seen on the scale of Figure 
5. The small protrusion on the drop, left by the breakup, quickly relaxes to an almost flat interface. 

The asymmetry of the breakup was already noticed in experiments |p5| , p6| . However, one must be careful not to 
apply our results to those experiments directly, since they are on length and time scales z' 3> l,|i'| ^ 1, far away 
from the asymptotic behavior. Still our similarity solutions could play a crucial role for the shape selection even in 
this "inviscid" regime, since all solution must ultimately match onto the asymptotic behavior. Clearly, an extension 
of our theory to the almost inviscid regime seems highly desirable. 
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Recently, an experimental study of drop formation in a highly viscous fluid has been reported | p7[ . Qualitatively, 
the shape of the interface adjoining the primary drop agrees well with Figure 5. Also, the length and time scales 
of the similarity solution, as given by the present theory, have been used to analyze the data and are found to be 
consistent with experiment. Unfortunately, at the times shown in Figure 2 A and Figure 2B of ||2^, the straining due 
to the falling drop is still appreciable compared with the scales of the similarity solution. Also, there is no independent 
measurement of t' available, which makes a meaningful comparison with theory difficult at present. We will discuss 
the process of repeated necking, reported in ||2^, below. 

Extensive experiments with high-speed jets, where gravity is irrelevant, are also in progress ||2^ . The stroboscopical 
method employed for example in ]29[ | allows to determine t' independently, so comparison with theory can be made 
without adjustable parameters. Preliminary results show nice agreement with theory before breakup. After breakup 
a quantitative comparison with theory is not yet possible, due to air drag on the rapidly receding neck, whose effect 
is not yet included in the equations. 

Therefore, we will use numerical simulations for a detailed comparison with theory. In particular, we would like to 
verify the prediction of the theory that the same similarity solutions are always approached, independent of boundary 
or initial conditions. Indeed, some simulations have already been performed on the breakup of a Navier-Stokes fluid 
pof , but they are not sufficiently close enough to the singularity to allow for a meaningful comparison. This is because 
in the asymptotic region Navier-Stokes computations become prohibitively expensive. To make simulations feasible, 
one has to resort to approximations. 

As model equations, we take the generalized form of the asymptotic equations (^0|) and (^). Extensive simulations 
of this system before breakup were already reported in |Q , , and . The equations read 

OtVo + vod^vo = d^p + g , (61) 



dtH + vod^H 



-{d,vo)H/2 



(62) 



where 



dlH 



H{i + {d,Hyy/^ (1 + (a,i/)2)3/2 



(63) 



So apart from the asymptotic terms already contained in ( p^ and ([ig|), (63) contains the exact expression for the mean 
curvature of a body of revolution. The system (|6l|)-(|63|) was supplemented with two types of boundary conditions 



In the "jet geometry" we fix the values of H and vq at two fixed points z+ and z_ 

H{z±,t) = H±(t) , 



(64) 



vo{z±,t) ^ v±{t) . (65) 

Hence here we envision a jet of length z+ — z- with nozzle radius — H- = ro and speed v+ = v- = V. At some 
point in time a small perturbation is applied to the speed V- at the nozzle and the jet breaks up according to the 
Rayleigh instability. The jet speed is so high that gravitational effects can be neglected, and thus g = 0- 

In the "drop geometry" fluid is released slowly from a tap. Thus at the opening of the tap,z_ say, boundary 



conditions (64) and (|65|) hold, while at the lower end of the drop the boundary moves with the fluid. This means we 
have 

H{z+{t),t)^0 (66) 

and 

vo{z+it),t)=dtz+{t) . (67) 

In this experimental situation gravity is of course important, as initially gravitational and surface tension forces are 
balanced, and the drop assumes an equilibrium shape ]3l| ]. These shapes are reproduced exactly by the stationary 
solutions of (|6^)-(|6^). Eventually, gravity overcomes surfaces tension and the drop falls and subsequently pinches off. 

The implementation of boundary conditions as well as the numerical procedure is explained in detail in . In 
[0 and 1^ simulations of (|6l|)-(l6^) have been used to reproduce experimental interface shapes both for high and 
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low viscosity fluids in different geometries. In particular in the case of a slowly dripping tap pq |, both boundary and 
initial conditions are known and comparison with experiment can be done without adjustable parameters. Thus the 



excellent agreement between simulation and the experimental shape of a falling drop at the pinch point |14 seems 
highly significant. Therefore we are confident that (^)-(|6^) represents a good approximation to the Navier-Stokes 
equation not only close to the pinch point, but also for earlier times and including the crossover to the boundary. 

We have performed systematic tests of the predictions of the present theory, in particular investigating the inde- 
pendence of the singular behavior near break- off from boundary conditions. For all runs, both in the jet and the drop 
geometry, and independent of the nozzle or tap diameter and of the viscosity, we always found the flow to converge 
onto the similarity solution predicted by the present theory. 

Figure 6 shows this convergence for a typical run in the jet geometry. The nozzle diameter is 100 in units of 
The solution near the singularity has been converted to similarity variables, thus giving (/)(^) and ^^{C) using the 
transformation (|2^). Shown is the predicted similarity solution as a solid line, and the computed solution at times 
\t'\ = 0.39, 0.13, 0.043, and 0.014, represented by dashed, chain-dashed, dot-dashed, and dotted lines. It can clearly 
be seen that the range of validity of the similarity solution expands like in the similarity variable ^. This 

means there is a fixed region in z' where the similarity theory applies, in agreement with the statements of Section 3. 
At the boundary of this region, the slope dz'h becomes of order unity, and the expansion in orders of e breaks down. 

The motion shown in Figure 6 occurs on scales widely separated from those imposed by the boundary conditions. 
The time distance from the onset of the linear instability to the singularity is tQ/tn — 32284, much larger than the 
relevant t' . Similarly, the nozzle diameter, converted to the similarity variable ^, is <^ = 157, 272, 473, and 829, for 
the times shown. Clearly the motion near the singularity has become independent of these imposed length and time 
scales. The same will happen for any boundary condition, as both the typical time and length scale shrinks to zero 
near the singularity. 

Next we test the convergence onto the similarity solution after breakup. Since there is a moving tip, we modify the 



viscous term in (61) to regularize the tip. The force balance now reads 



dtVo+VQdzVQ^ ~-dzP + iv — ^ - g , (68) 

p 

where a is a free constant. By varying a we can test our prediction that the shape of the interface after breakup does 
not depend on the regularization employed. 

To produce an initial condition after breakup, we take a simulation before breakup, which has progressed to a 
time distance of \t'\ = 10~* from the singularity. Then we cut the solution at the minimum and interpolate H to 
zero with a polynomial, so as to keep the highest derivatives smooth. This we take as the new initial condition after 



breakup and let the solution evolve under (68), (|62|), and (|63|). For a wide range of values of a in (|68|), we always find 
the solution to converge onto the similarity form found in Section 4. Figure 7 illustrates this convergence for a run 
which has the same boundary conditions and material parameters as the one shown in Figure 6 before breakup. The 
constant a was chosen to be 1. Again, solutions were converted to similarity variables. The full line represents the 
predicted similarity solution, the dot-dashed and the dashed lines show the numerical simulations for \t'\ = 0.006 and 
\t'\ = 0.06 after the singularity. The dotted line is the similarity solution before breakup, shown for comparison. It 
can clearly be seen that after the solution has been cut in two halves it rapidly converges onto the predicted similarity 
form. This is independent of both the regularizing term in (|6^) and the procedure by which the solution is cut. 

Hence both before and after the singularity, we have always observed convergence onto the similarity solutions if |<'| 
is small. Still it would be very useful to have a better mathematical understanding of the approach of the similarity 
solution for the full Navier-Stokes dynamics. Even for the simplified model equations (|6l|)-(|6^) the convergence we 
found numerically is far from being a trivial result, as there are higher order derivative terms like d'^H coming from 
the pressure. In principle, although these terms are multiplied by a small number t' close to the singularity, they could 
make a singular perturbation, which changes the asymptotics. However, it is well beyond the scope of this paper to 
explore these questions in detail, so at present we have to rely on the ample numerical evidence. 

Another important question is the stability of the similarity solution to small perturbations. This has been studied 
in the framework of the asymptotic equations (p^,(p^) in |^2|, both numerically and analytically. The result is that 
the similarity solutions are linearly stable as expected, since there are observed numerically. On the other hand, they 
are unstable to finite amplitude perturbations of wavelength comparable to the minimum radius of the fluid neck. 

As soon as a perturbation is large enough, it will start to grow and eventually forms a new similarity solution with 
its own zo and to. For any finite number of such perturbations, the singularities are separated in space and the present 
theory strictly applies. However, if one explicitly adds an external white noise source to the Navier-Stokes equation, 
perturbations are introduced on all time scales arbitrarily close to the singularity. This allows for the appearance of 
a "rough" interface as described in and |^^, consisting of an infinity of interacting similarity solutions. Locally, 
the form of each of those solutions, seen as "necks" in experiment, is consistent with the present theory. 
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In ||3^, a threshold length scale 



-Cf-^j (69) 

was identified, below which thermal fluctuations become important. The relevant thermal length scale for surface 
perturbations is It = {ksT /^Y^'^ . For a mixture of 85% glycerol and 15% water ithres is l/im. Thus, in the presence 
of thermal fluctuations, the microscopic length scale imicro introduced in Section 3 may be replaced by ithres- for Hmm 
larger than ithres the Navier-Stokes equation is applicable, on smaller scales the equations are inherently stochastic. 

An obvious benefit one expects from the universality found in the present paper is the unique continuation of 
Navier-Stokes simulations through the singularity. One slight problem lies in the nonanalyticity of (f)^ and ip^ if one 
wants to use similarity solutions as new initial conditions after breakup. Although the pressure itself would be finite, 
pressure gradients and the viscous term would be infinite at the tip. Thus it is better to use regularized similarity 
functions, where a small but finite cut-off parameter £ has been introduced. The resulting initial conditions for the 
new Navier-Stokes problem after the singularity would be arbitrarily close to the similarity form, but still finite at 
the tip. 

In conclusion, we have shown that the Navier-Stokes equation carries us through the bifurcation point where at first 
it seems meaningless. As usual, classical hydrodynamic theory has a much wider range of applicability than purely 
microscopic considerations would tell us. 
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a+ 
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K 


-1.5699 


0.030432 


-1.6024 


0.030426 


-3.066 


4.635 


0.0723 


6.047 X 10~* 


57.043 



TABLE I. Some characteristics of the similarity functions (j)'^ ,tp^ before breakup. The symbols and 00 stand for the 
position of the stagnation point, where the fluid is at rest in the frame of reference of the interface, and the radius of the 
interface at that point. The minimum value of is 0min, and ^min is its positiou. The function ip'^ reaches a maximum value 
of ipmax- The numbers and stand for the limits limj-^ioo 0^(C)/C^ and lim^^±oo ^^(C)Ci respectively. All numbers are 
accurate to the decimal places shown. 



^neck 




^drop 




17.452 


0.06183 


0.4476 


0.6180 



TABLE II. Characteristics of the similarity functions ,tp~ after breakup. The tip position of the left hand, or neck side 
is Cnecfe, and the expansion coefficient 01, cf. (^^, is 0necfc. Correspondingly, ^drop and (pdrop uniquely determine the "drop" 
side of (p~ and The values of and are the same as before breakup, cf. Table 1. 
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FIG. 1. A sketch of the flow geometry investigated in the present paper. The radius or "height" of the free surface at a 
point z on the ajcis of symmetry is Hiz). The velocity fleld inside the fluid is v(2;,r) = Vz{z,r)&z + i'r(z, 7')er. 



FIG. 2. A plot of the similarity functions (^"'", (a), and ■i/'^, (b), before breakup. Note the strong asymmetry. 



FIG. 3. A cartoon of a receding tip after breakup. The position of the tip is z'^^pit). 



FIG. 4. The similarity functions 0", (a), and t/j", (b), which are unique continuations of (^i^ and i^)'^ to times greater than 
io. The asymptotic behavior for ^ — > ±oo is by definition the same as before breakup. On the left is the rapidly receding 
"neck" part of the solution, on the other side is the drop. The points at ^necfe and ^drop, from where the interface is plane, are 
marked by diamonds. 

FIG. 5. The breakup of a mixture of 5 parts of glycerol in 4 parts of ethanol, as calculated from the similarity solutions. 
Part (a) shows three profiles before breakup, in time distances of 46/xs, corresponding to \t'\ = 1,0.55, and 0.1. In part (b) the 
same is shown after breakup. 

FIG. 6. Simulation of (pl|)-(]6^ in the jet geometry. The profiles close to pinch-off were converted to similarity variables. 
The full line is the prediction of the present theory; the dashed, chain-dashed, dot-dashed, and dotted lines represent the 
simulation at \t'\ = 0.39, 0.13, 0.043, and 0.014. The inset contains a blowup of the central region with only the latest time, 
\t'\ = 0.014. 



FIG. 7. The approach of the similarity function 0" by the solution of (|68|),(|62|) and (|63|) in the jet geometry, transformed 
to similarity variables. The fiuid neck is severed at \t'\ = 10~* before breakup. The full line is , the dotted line the solution 
before breakup. The dot-dashed and the dashed lines show the simulation at \t'\ — 0.006 and \t'\ = 0.06, respectively. 
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Figure 2b: 
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Figure 4b: 
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Figure 5a: 
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Figure 6b: 
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